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Abstract 

This paper describes an alternative method of generating fixed points of 
certain substitution systems. This method centres on taking infinite words 
consisting of one repeated letter per word. These infinite words are then inter- 
laced to form a new, more complex, infinite word. By considering particular 
limits of interlacings of words, fixed points of substitutions are generated. This 
method is then extended to two dimensions, where a structure equivalent to 
a well known aperiodic tiling (the Robinson tiling) is constructed. 



Introduction 

In the first section of this paper we describe our fundamental operation, the D 
operation. This operation acts on two infinite words (or ID arrays) to produce a 
new infinite word, by placing a letter of the second word between every two letters 
of the first word. We then show how this operation can be linked to Toeplitz 
sequences, and show how iterating the operation can create more complex words. 
In the second section of this paper we introduce a generalisation of the fl operation, 
namely fl n . This operation places a letter from the second word after the (kn) th 
letter of the first word (for all fceZ). We then consider limits of these operations, 
and prove that carefully selected limits not only exist, but are fixed points of specific 
substitutions. Hence we have an alternate method of calculating fixed points of (a 
family of) substitutions on words. The results in the first two sections link in to 
known Toeplitz sequence theory. The third section extends these results to two 
dimensions, detailing how this method can generate a Robinson tiling (one of the 
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first 'aperiodic' tilings, as described in [S]). The paper ends with a brief summary 
of possible uses and avenues for further investigation. 

This paper attempts to err on the side of clarity, instead of brevity. As such 
we have included several worked examples were we walk the reader through new 
concepts slowly. Of course the reader is free to skip worked examples. 

On a first reading we suggest that the reader briefly skims subsections 1.1 and 1.2 
and the start of section 2, then skips ahead to subsection 2.4 in order to understand 
what the method can do in the one dimensional case. After this, section 3 can be 
read to see how the method extends to two dimensions, or the skipped subsections 
can be read for concrete proofs. 

1 The n operation 

The R operation uses the concept of an Z° '-array. 

We will use these arrays to describe both the 1 and 2 dimensional cases, using 
similar terminology for both dimensions. In this section we will introduce the (1 
dimensional) fl operation, leaving the 2 dimensional case to the third section. Let 
us consider the following definition of Z d -arrays. This formulation of Z d -arrays was 
first studied in j2], and extended in [1]. 

Definition 1 (Z d -array). A Z < -array (with alphabet S) is a function 
Z : Z d ^ £ 

The precise formulation of Z d -arrays we are using is new to this thesis, but is 
equivalent to the version found in [TJ. 

Definition 2. Choose a set of basis vectors {v±, . . . , vj} of M. d . Define the underlying 
lattice Latt as {a\V\ + 02^2 + • • • + adVd\a\, 02, . . . G Z}. Define p : Z d H- Latt as 
p(ai, 0,2, ■ ■ ■ , ad) = a\V\ + 02^2 + • • • + adVd- This is clearly a bijective function. 

Then the array on a lattice (Z, Latt) is the function; 

Z 1 : Latt 1 — y E such that Z 1 o p = Z. 

In other words, Z is the pullback of Z 1 . Note that if you have an array Z and 
an underlying lattice Latt, you can calculate Z'. Similarly, if you have Z ! you can 
calculate Z. Thus we will often refer to 'arrays on a lattice' as 'arrays', when the 
lattice is fixed. 
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These definitions correspond to assigning an element from an alphabet E to 
every point on a lattice embedded in M. d . The underlying lattice Latt encodes which 
lattices we choose and the Z d -array encodes what elements from the alphabet get 
assigned to the particular lattice points. From now on we will use the convention 
that when no underlying lattice is mentioned, we are using the integer points in M. d 
for our underlying lattice. In one dimension this corresponds to bi-infinite sequences. 

By considering the Voronoi diagram of a lattice, we can construct an associated 
tiling for any lattice. A tile in this case would be the Voronoi cell of a point p in 
a lattice A C M d (namely the set of all points x G W 1 such that d(x,p) = d(x,A)). 
Thus we can consider the D operation as an operation on tilings as well. 

1.1 1 dimension 

The major advantage of studying the ID case is that we can express the array as a 
bi-infinite sequence, enabling us to use results applicable to sequences. Thus we will 
define the D operation for ID arrays, and Toeplitz sequences, a class of examples 
originating in the field of dynamical systems [5] [?]. 

Definition 3 (An B, the Superposition operation). Let A and B be Z-arrays, 
namely A : Z ^ E, B : Z ^ E. 

Then define A D B as follows; 



Definition 4. Define fx(A) as the function sending an array A to Af] X, where X 
is another array. 

1.2 Motivation 

We can consider A D B in one dimension as the overlaying, or merger of two separate 
tilings. Imagine the tilings A and B as two infinite translucent physical strips with 
the strip representing each tiling. The operation A D B corresponds with placing 
A down over the origin normally, then shifting it half a unit to the left. You then 
place the strip for B down over the origin normally, effectively interlacing points 
from A and B. Looking through A and -B's strips, the tiling A D B can be read 
off. Note that the distance between consecutive points needs to be scaled back up 
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to 1 (by expanding about the origin by a factor of 2). See picture [T] for a schematic 
motivation of how to create A (1 B from two tilings A and B. 




Figure 1: A schematic motivation for Af]B. 

We will now give the definition of a Toeplitz sequence. We will also define a 
'null' array, a concept used in the field of Toeplitz sequences. Note that there are 
many equivalent definitions of a Toeplitz sequence. We will use one more easily 
applicable to this thesis, from the reference [5] (with a very similar format to [1]). 

Definition 5. Let A be a finite set of at least two elements. Let A* be the set of 
finite sequences, or words over A. If w G A*, let \w\ denote its length. Let Q = A z . 
If S G f2, n G Z and p > 1, then let 

S n S n +l ■ ■ ■ »5n+p-l 

denote the word of length p appearing in S starting at position n. Thus S n is the 
nth letter in the sequence. 

Definition 6. An element S G f2 is called a periodic sequence with period p G N if 
S t = S t+P for all t G Z. 
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Definition 7. An element S G f2 is called a Toeplitz sequence if it is not a periodic 
sequence, and satisfies the following condition; 

Wn G Z, 3p > 2 such that V/c G Z, S n+ k p = S n 

Definition 8. A p-periodic part of a bi- infinite sequence S G A z is; 
Per p (S) = {n G Z : Vfc G Z, = S n } 

Remark. Consider a sequence T, which is not periodic. If every point T fc G T is in 
a p t -periodic part for some p t , T is a Toeplitz sequence. This is because if a point 
T n belongs to a ^-periodic part, then by definition of periodic part, T n+ k p = T n . 

Definition 9. An almost Toeplitz sequence is a sequence where all but a finite 
number of points are in a periodic part. 

In this chapter we will be creating (almost) Toeplitz sequences. The points not 
in a periodic part will be a finite word directly to the right of the origin. We will 
refer to this word as the seed, for reasons which will be apparent later. 

Definition 10. A null array * is the unique array with alphabet {*}. 

Regarding the D operation, let us define it to be an operation that is conven- 
tionally evaluated, as follows; 

Definition 11. A n B n C := ((A nB)nC) Furthermore; A x n A 2 H . . . n A n+l : = 

(Ai n A 2 n . . . n A n ) n A n+i ) for all neN. 

Definition 12. (A fl B) 2 = An B n An B 

if a o f B ) 2 (x) = f A °f B °fA° f B (x) = xnBnAnBnA 

Definition 13. If Bi, B 2 , B s , ... is a sequence of arrays Z H- S, say they have the 
limit B : Z H- S if Vn G N 3m such that = Bi(x) for all — rt < x < n and 

i > m. 

For studying infinite fl operations, we need to be careful when looking at the 
limit. 

Consider two arrays, and 1, where 0(f) = and l(v) = 1 for all v G Z. Now 
consider the series {X i \ c *l 1 — 0, D 1, D 1 fl 0, . . . as shown in figure El 

By examining definition [3], the reader can observe that points in the even posi- 
tions of an array An B are only dependent on the array B. Thus the same result 
applies for a function fg(A), since it is merely different terminology. 
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Figure 2: The first few points of 0, D 1 and other arrays. 
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We will use this fact to show that no limit of the series {X i }°? =1 exists. 
The sequence {X i \ c *l 1 = 0, D 1, D 1 D 0, . . . can be rewritten as 

0, /i(Q), /o(/i(Q)), . . . , / a (X 2i+1 ), fi(X 2(i+1) ), fo(X 2{l+1)+1) ) . . . 

More precisely, X 2i = fi{X 2 i-i), X 2i+1 = fo(X 2i ) and X 1 = 0. 

The value of even positions in X 2 i = fi(X 2 i~i) is 1. The value of even positions 
in X 2i+ i is 0. Thus the value of points in even positions will be different in X 2i and 
X 2i+ i, for any % 6 N. Thus there does not exist a limit of this sequence. 

Limits of D operations do not (usually) exist, but we can instead use the com- 
mon idea of convergent subsequences, in an attempt to get interesting results. For 
example with the limit ofOninOnin... we could consider the limit of 

oni, (on i) 2 ,(on i) 3 ,... 

or 

o,onino, (oni) 2 no, ... 

(These subsequences do have properly defined limits, which will be proven in a 
later section). Considering different subsequences may give us different values for 
a limit. Note that the work in this section implies that only superpositions of the 
form A\ n A 2 n . . . n A n n A n A n A . . . could have limits when considered in the 
naive way. 



2 An n B 

We will now define an operation without obvious equivalent in higher dimensions, 
n n . We will then describe the limit of a sequence of n n superposition operations. 

Definition 14. Let A and B be one dimensional sequences. Define An n B : Z H- S 
as follows; 

(An. B|W J fl ^ if „ = *(,, + 1) fa tez, 

I A[v — L^TiJ) otherwise. 

Intuitively you place Bq G B at the origin, bumping Aq left one unit, and this 
bumping all points left of the origin one unit left. You then add £>i to the right of 
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An, bumping tiles A^ for N > n right by 1. Similarly you add B_i to the right of 
A_ n , bumping tiles away from the origin to make space. You then place B 2 to the 
right of A2n and B^2 to the right of A_2 n , and iterate for B3 and -B-3, and so on. 
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Figure 3: The arrays A, B, and A ff B for 1 < i < 5 

Definition 15. Define fx{A) as the function sending an array A to an array An n X. 

Remark. Note that the fl 1 operation is equivalent to the D operation defined pre- 
viously. 

2.1 (runs) 00 

We will now consider a subset of infinite compositions of the fl and fl n operations. 
Note that the fl operation is equivalent to the operation fl 1 . Thus, in the 1 dimen- 
sional case, fl n can be considered to be a more general operation. We will start with 
the basic definitions, before giving illuminating examples. 

Definition 16 ((f] s A fl* B) m ). Consider two sequences A and B. Let s,t,m E N. 
Then define (f] s A n* B) m as follows; 

(nM n* B) m ■= A n f b n s A n f b n s nMn'B 

2m. 

Similarly (for sequences Ai, and integers EN). 
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(n ai Ain a2 A 2 . . .n an A n ) m ■= A± n a2 A 2 ... n a - A n n ai A n a * A 2 rT"- 1 A n _ x n a - K 

mn 

Definition 17 ((fl s An' B)°°). Let A and B be sequences. Let s,t G N. Then 
define; 

(n s A n* B)°° := lim {n s A n* B) m 

m— >oo 

if the limit exists. 

Similarly (for sequences A iy and integers Oj 6 N), 

(n ai A! n a2 A 2 . . . n a - ^ n )°° := lim (n a M! n a2 A 2 . . . n a - A n ) m 

if the limit exists. 

We will now work towards calculating the limit of repetitively applying a series 
of (l n N operators, with different periodic arrays N. We will eventually prove that 
you can create substitution sequences by this method. To improve clarity, we will 
run through a simple example (with associated proofs) first, to give the reader some 
intuition into this 

2.2 0D1 and its limit (0 n 1)°° 

Recall that (0 n l) m is defined as 0nin0ni...n0ni . Similarly, (0 n 1)°° is 

2m 

defined as the limit of (0 D l) m as m — > oo. 

We will show that (0 Dl) 00 exists, and is equivalent to an (aperiodic primitive) 
substitution tiling. 

Consider the sequence of arrays *, /i(*), /o ° ■■■,(/o° fi) m (*),fi ° (/o ° 

fi) m (*), . . .. Part of this sequence is shown in figure HI Recall that * is the null 
array. 

As we repetitively apply the fo ° fi operation to our array it appears that larger 
and larger patches of the array are independent of what the starting array was. We 
will now make this rigorous. 
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Figure 4: Sequence of arrays. 

Definition 18. Consider an array Ai fT 2 A 2 fT 3 A 3 . . . fT m A m . A point p in the 
underlying lattice of A 1 fT 2 ^4 2 • • • H'™" 1 A m is called an undefined point if * fT 2 ^4 2 H'™ 3 
A3 . . . fl n " 1 A m maps that point to *. 

Thus from definition [3], any array of the form X fl 1 (ie fi(X)), where X is any 
array, has undefined points at odd positions. 

Lemma 1. Let m G No and k G N. If a point is defined in f\o (/o /i) m (*), it will 
be defined, and take the same value, in fi o (/q fi) m+k (*)- 

Proof. Consider /o(^0; for any array X. By the definition of the fl operation, we 
know that fo(X) will have defined points at positions 2s (s G Z), with value 0. 
(Similarly, fi(X) will have defined points with value 1 in the same positions). 

Consider f^of^X) which is equivalent to fi applied to fo(X). By the definition 
of the fl operation, a point at position 4s — 1 in fi o fo(X) will share the same 
value as a point at position 2s in fo(X). Thus j\ o fo(X) will have defined points 
at positions 4s — 1 with value 0, and defined points at positions 2s with value 1. 

Let us generalise to an array o f i2 o . . . o f in (X). This array will have defined 
points at positions 2 n s — 2 n_1 + 1, with value i n . We will use induction on n to prove 



10 



this. 

Regarding the initial case, fi^X) has defined points at positions 2s as required 
(via definition [3]). For the induction step, consider f i2 of i3 o...o f in+1 (X). This 
will have defined points at positions 2 n s — 2 n ~ l + 1, with value i n +\, from our 
assumption. Call this array X n . From definition [31 fi 1 (X n ) will have a point with 
value X n {y) at position 2v — 1. Thus the set of points with value i n+ i in fi^Xn) 
are {x\x G 2(2 n s — 2™" 1 + 1) — 1, s G Z}. This can be rewritten as S n+ i = {x\x G 
2 n+1 s — 2 n + 1, s G Z}, as required. Thus our induction holds. 

Thus /io(/o°/i) m+fc (*) will define every point /io(/o°/i) m (*) defines. This is be- 
cause /i°(/g°/i) m (*) will have defined points (with associated values) in |Ji=i m+1 &i 
and fi o (/o ° fi) m+k (*) will have defined points in (j*-2( m + fc )+ 1 Furthermore, the 
defined points will have the same value, since the sets S n = {2 n s — 2 n ~ 1 + l|s G Z} 
are disjoint for n > 1, thus the value of a point is solely determined by which set Si 
it is contained within. 

For an explicit proof that the sets S n are disjoint, consider S a and Sb, for a, b G N, 
a b. Assume WLOG that a < b. Then b = a + c, for c > 0. The sets 
and Sb will have non-zero intersection if there exists an integer solution to the 
equation 2 a s - 2 a ~ l + l = 2 h t- 2 6 " 1 + 1. This can be rewritten as 2 a s - 2 a ~ 1 + 1 = 
2 a+c t _ 2 a+c-i + 1; and sim piifi e( i to 2s - 1 = 2 c+1 t - 2 C . The left hand side of this 
equation can only take odd values, and the right hand side can only take even values 
(since c > 0). Thus S a and Sb are disjoint. 

□ 

Theorem 1. lim„ woo (0 D l) m exists. 

Proof. Note that another equivalent way of writing lim m _ s>oo (0ni) m is as linv^oof/io 
(/ooA)" 1 " 1 ®]. 

Consider the lattice point at position '1' in any array. For any array X, the fx 
operation does not change what value is assigned to this point. An array (0 D l) m 
is constructed by taking a array, and applying multiple f\ and fo operations to it 
in turn. Thus the lattice point at position '1' in an array (0 D l) m will be a point of 
type 0, for all m. Thus the lattice point at position '1' is well-defined in the limit. 

For our next step, we will show that in the limit, the only undefined point is 
at position '1'. Call this point the seed. Consider /i(*). The points at positions 2 
and are defined. On applying fg to this array, the point at position 2 is shifted to 
position 3, and the point at position is shifted to position —1. Thus in fgo /i(*), 
the points two units away from the seed point are defined from fi, and the points 
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one unit away from the seed are now defined by Jo- Every time a new fa function 
is applied to the array, the patch of defined points grows by (at least) one. Thus in 
the limit, every point (barring the seed point) is defined for some f\o (/o o fi) m (0). 
From the previous lemma, we know that every defined point will keep the same 
value as m —> oo. Thus these lattice points are well-defined in the limit. Thus 
lim m _ > . oo (0 fl l) m exists. 

□ 

Theorem 2. (Ofl 1)°° is an almost Toeplitz sequence, with the only undefined point 
being at position '1'. 

Proof. We will use induction to show that every defined point in (0 D 1)°° is in 
some periodic part. Construct the sequence {Xi} where X\ = /i(0), X 2 = fo(Xi), 
X 2i = fo{X 2i -i) and X 2i+1 = fi(X 2i ) for i > 2. 

Consider the defined points of X\ = /i(0). These points occur at even positions, 
and thus form a 2-periodic part, as defined in definition |HJ Thus every defined point 
in X\ is in some periodic part. (In fact, every point is in a periodic part). 

Assume every defined point in X n is in some periodic part. Consider X n+ \ = 
fi(X n ) (where /j is either /q or fi, based on parity of n). (Denote the points of X n 
as x' t , and the points of X n+ i as x t , for t 6 Z). 

Any defined point in fi(X n ) will either correspond to a defined point in X n , or 
will be an even point. If it is an even point, it will be assigned a value of 1 (if n is 
even) or (if n is odd). Since all even points are assigned the same value, the even 
point of fi(X n ) form a 2-periodic part. 

Any defined point x 2t -\ G fi(X n ) will correspond to a point x' t in X n . By 
our assumption, the point x' t in X n belongs to some p-periodic part, for some p. 
This periodic part consists of a series of points x' t+kp with the same label (for k G 
Z). Therefore the p-periodic part will be mapped to a 2p-periodic part in fi(X n ), 
consisting of a series of points x 2 ( t +kp)-i with the same label. Thus any defined 
point in A n+1 belongs to some periodic part. By induction, any defined point in 
Xi, for any i, must be in a periodic part. We know from the last theorem that if a 
point is defined in Xi, it will have the same value in all sequences X i+2 k. Therefore 
once a point is defined (and thus belongs to a periodic part), it will always belong 
to a periodic part. We can therefore conclude that every defined point in the limit 
(0 H U°° i s m a periodic part. 

□ 
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Remark. Any finite composition of periodic arrays under the PI operation will be 
periodic. In particular, the sequence (0 D l) m is periodic, with period 2 2m ~ 1 . 

As a sketch proof, let {Xk}^L be any arrays with only one letter in each alpha- 
bet, possibly different for each array. (A similar proof applies for periodic arrays, but 
requires more complex notation for little added clarity). Consider some sequence 
F\ = /xi(Ao), Fi = fxi(Fi-i). Consider a point in X . It will be in a 1-periodic 
part. Under the fx x operation, this 1-periodic part will be mapped to a 2-periodic 
part. Any point in F\ will either be undefined (and hence in the 2-periodic part 
inherited from X ), or it will be defined. If it is defined, it will be in another 2- 
periodic part, this one consisting of points from X\ in even positions in the array. 
Thus F\ has period 2. 

Applying the fx 2 operation to Fi will send the 2-periodic part representing 
undefined points in F\ to a 4-periodic part in F 2 . The 2-periodic part representing 
points first defined by fx 1 will be sent to a 4-periodic part, and a new 2-periodic 
part will be created for the points from X 2 which have been placed in even positions 
in the array F 2 . Thus F 2 will have period 4 (the 1cm of all the periodic parts). 

Similarly F3 will consist of two 8-periodic parts, a 4-periodic part and a new 2- 
periodic part containing values from X3. In general, Fi will contain two 2*-periodic 
parts, and one 2- ? -periodic part for each < j < i. Thus Fj will have period 2\ 
(0 fl l) m can be expressed as F 2m _i for X 2 k = 0, X 2 k+i = 1. Thus the result follows. 

Theorem 3. The sequence (0 (~l 1_)°° is a fixed point of the substitution, 

a :0 H> 0101 
1 ^ 1101 

More precisely ff((0n If ) = (On If . 

Furthermore lim m _ 5 . oo ((0_nl) m n0) is also a fixed point of a substitution, namely; 

a:04 0010 
1 ^ 1010 

More precisely a(lim m ^ oo ((0 n l) m n0)) = lim m ^ oo ((0 n l) m n 0). 
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Proof. First consider lim((0 D l) m HO). Note that this limit exists, since it is equiv- 
alent to lim/o((On l) m ), and lim((0 n l) m ) is well defined. 

We know that lim m ^ oo ((0 (~)l) m DO) must be invariant under fo o /i, since fo° fx 
is equivalent to increasing the value of m by one. From theorem [2} we know that 
the point at position '1' of lim((0 D l) m DO) must be of value 0, since neither fo or 
fx can alter the value of the point at position '1 ', and it is of value '0' in the array 
(the first sequence we are building the limit from). Denote the point at position 
'1' the seed point. 

Examining the proof of theorem [21 note that points defined by the first ap- 
plication of fo must be in a 2-periodic part, and points defined by the following 
application of fi will be in a (2 x 2)-periodic part. Thus we can deduce that the 
highest periodic part in fo o is of period 4, thus fo o is periodic, with 

period 4. Furthermore the repeating period is the word * 1 0. 

Now consider a generic sequence t = . . . totfa • ■ ■> and what points are defined 
from it under the application of fo ° fx- The seed point, ti gets sent to the word 
tx 1 0. Via the periodicity of fo ° we know that the point t 2 is sent to the 
word t 2 1 0, which is joined on to tx's word by concatenation. Similarly, t 3 's word 
is concatenated to the end of t 2 , s word, and so on for all £j. This is of course the 
function a from the definition of the substitution. Since the seed point is fixed for 
all m > 1, the points generated from it are fixed for all m > 2. Similarly, all points 
generated from those points are fixed for m > 3, and so on. Thus all points defined 
via the application of fo o fi are fixed in the limit. 

Since every point in the underlying lattice is defined by fo°fx(t), we can conclude 
that lim((0 n l_) m D 0) can be generated via a substitution. Since fo o /i(lim((0 D 
l) m n 0)) = lim((0 n l) m n 0), it is a fixed point of that substitution. 

For (On!) 00 use a similar proof, using the fact that (On!) 00 is invariant under 
/Wo- □ 

We will now generalize these proofs to arrays of the form (n ai Aid 0-2 ^n" 3 . . .n a " 
An) 00 , where A\ are periodic arrays. 

We will also show these types of arrays are (almost) Toeplitz arrays. 

2.3 (n fli ^i rY 2 A 2 n° 3 . . . n a " A n )°° 

We can extend the proofs of lemma [1] and theorem [1] to cover the more general case, 
(n ai ^4i H a2 A 2 n a3 . . . n an A n )°°. To avoid excessive duplication, we will sketch how 
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to alter the proof to the more general case. 

Theorem 4. (n ai Ai (T 2 A 2 H as . . . H a ' 1 A n )°° is well defined. 

Sketch proof. Take a = min{ai, . . . , oat}. Any point in an array X with position 
strictly between and cr is invariant under any operation n ai Ai. Thus the first 
(cr — 1) points of Ai form a 'seed' which is invariant, and thus the values of the seed 
points are well defined in the limit. 

Ignoring the seed, the first undefined point in fl ai A\ fl a2 A 2 n as . . . fl an A n cannot 
be closer to the origin that position (cr + 1), since the first (cr — 1) points are in 
the seed, and a new value will be inserted at position cr during every iteration of 
n ai Ai (l a2 A 2 (l a3 ... (l an A n , by the definition of a. 

Similarly, the first undefined point of (n ai Ai fT 2 A 2 (T 3 . . . H a " A n ) 2 cannot be 
closer to the origin than cr + 2, and in general the first undefined point of (n ai Ai fl a2 
A 2 fl a3 . . . rT n A n ) m cannot be closer than (cr + m). 

(n ai Ax n a2 A 2 n a3 ... (l an A n ) m can be rewritten as; 

(fZ ° fZ~-\ ° • • • ° fZ) k (( naiA i • • • na " A nT~ k 

for m > k. 

Thus any defined point in (n ai Ai fl a2 A 2 n as . . . fT" A n ) k will have the same value 
for any (n ai An a2 A 2 n a3 . • n an A n ) m for all m > k. Thus ((T 1 A 1 n cl2 A 2 n a3 . . .n an A n )°° 
is well defined. 

□ 

Theorem 5. If A { are periodic arrays, then (n" 1 ^ fl a2 A 2 fl a3 . . . fl a ™ A n )°° is an 
(almost) Toeplitz array. 

Proof. Note that if Ai is a periodic array, every point in A, is in a periodic part. 
The maximum period of these periodic parts is equal to the period of the array Ai, 
denoted p(Ai). 

We wish to show that every defined point in (n ai A 1 n a2 A 2 n a3 . . . H a " A n ) m is 
in a periodic part (ignoring the seed). Construct a sequence C = {0(2)}°^ of arrays 
where the sth term in the sequence consists of the array A\ with s f^. operations 
applied to it in ascending cyclic sequence. 

The first few terms of this sequence are; 

c(l) = A x n a2 A 2 . 
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c(2) = A t n a2 A 2 n a3 A 3 

c(3) = A 1 n a2 A 2 n a3 A 3 n a4 A 4 

c( n - i) = A x n a2 . . . n a " A n 

c(n) = Ai n a2 . . . n a ™ An n ai Ai 

c (n + i) = Ai n a2 . . . n a " A n n ai A 1 n a2 A 2 



Explicitly c(s) is; 

. . .where j is such that s = (n — 1) + (k — l)n + j, and j < n. 

The arrays (fl ai Ai fl a2 A 2 fl a3 . . . fl a " A n ) m form a subsequence of this sequence. 
Thus if we can show that every defined point in any array (from C) is in a periodic 
part, ignoring the seed, then we are done. We will use induction, showing that this 
statement holds true for c(l), and that if the statement is true for c(s) it is true for 
c(s + l). 

For the initial step of our induction, consider /^ 2 (*), alternatively known as 
* fl a2 A 2 . The points defined by belong to a finite number of periodic parts, with 
maximum period (a 2 + l).p(A 2 ). Thus every defined point in f^(Ai) belongs to a 
periodic part. 

For our induction step, assume every defined point in the array c(s) is in a 
periodic part. Without loss of generality, take any one of these periodic parts, p, 
of period x. Consider Ja^- Consider a section of c(s) of length lcm(x, Oj+i). The 
operation f A 3+ \ adds a point every aj + i units. Thus under the f/ + \ operation any 
section sec of length lcm(x, ctj+i) is expanded to a section sec of length 

. . lcm(x, a,+i) 

/ = lcm(:r, a j+1 ) H v 3 ' . 

a j+1 

Furthermore, after the section of length lcm(x, a J+ i) has been expanded, the 
following section in c(s) of length \cm(x, a,j + i) will also be expanded to a section 
of length /. Via the properties of the 1cm function, the positions of points from p 
will be the same in sec, the following section of length /, and all following sections 
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after that. Thus a x-periodic part in c(s) will be turned into a /-periodic part in 

Thus by induction, every defined point (not in the seed) in (n ai Ai fl a2 A 2 fl a3 
. . . n a " A n )°° is in a periodic part. By theorem OJ (C\ ai A 1 n a2 A 2 H as . . . D a " Ai)°° 
is (almost) Toeplitz. □ 

We will now introduce a natural extension of the well-known concept of a sub- 
stitution. The major change will be that the defining function of the substitution is 
based on words, not letters. 

Definition 19. Let A be an alphabet, A* be the set of finite words over A, and 
W C A*. Let A°° be the set of bi-infinite words, and A°°w be the set of bi-infinite 
words with unique decomposition into a sequence of elements from W. Let A*w be 
the set of finite words with unique decomposition into a sequence of elements from 
W. 

The function a : W (-> A* is a valid substitution on words if it satisfies the 
following two properties; Firstly it must induce a well defined function a°° : A*w 
A*w. Secondly for all w G W, there exists N such that \w\ < \<7 n (w)\. 

A cr-substitution tiling is an element A G A°°w when A C Im{a (X ') N for all 
N G N. 

For example a function a defined as cr(01) = 0110 , c(10) = 1001 would induce 
a function a°° sending 011010 ... to 011010011001 . . .. 

Theorem 6. (n ai A 1 n a2 A 2 n" 3 . . . n a?l A n )°° is a fixed point of a substitution on 
words. 

Proof. For brevity, we will define the function o f^~_ L o . . . o as F. Consider 
F(*). By the proof of theorem points defined by f%. will be in a pj-periodic part, 
where pi is some finite number. Thus, since we have only applied a finite number 
of f£. operations, there are a finite number of periodic parts. Thus F(*) is periodic 
(with period less than or equal to p\ ■ p 2 ■ ■ ■ ■ ■ p n )- 

Since F(*) is periodic, there must be a repeating sequence S of minimum length 
(starting at the origin). For example, let K{ be an array with alphabet {ki}. Then 
fx 2 ° w °uld have a repeating sequence * k 2 * k 2 ki k 2 . 

By theorem HI there exists a finite number of 'seed points' near the origin that 
are unchanged by f^., for any i G N. Call the word formed by these points w = 
WiW 2 . . . Wk- Consider the number n of undefined points in S. If \w\ > n, then we 
can repetitively apply our function F to the first n points of w to define the entire 
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tiling. This is equivalent to a substitution rule sending wi . . . w n to the first \S\ terms 

of F(w) where w = W\W 2 ■ ■ ■ w n * * * For our earlier example of fx 2 o /#(*), 

this would give a substitution rule ui\w 2 w\ k 2 w 2 k 2 k\ k 2 . 

If \w\ < n then repetitively apply F to w until the first n points are defined. 
Then apply F to these n points as before. 

□ 

Let us consider some examples, in order to improve our intuition regarding this 
process. 

2.4 Examples 

Define X_ as the array where all points in the array have type X, where X can be 
any fixed symbol. 

Example 1 ((DO n 1 n 2)°°). 

(/a fi /o)(*) = h * * . . .) 

= /a(*101 * 101 . . .) 

= *2120212* 2120212... 

The repeating sequence is thus *2120212. Thus the substitution is x\ 1— >■ Xj2120212. 

Tne /irsi array m (DO Dl n2)°° is the array, and only the point at the origin in 
is unaffected by fo, fx and f 2 . Thus the seed of (flOnl n2)°° is the word ' ' ' 
has one point in it, which is equal to the number of unknown points in *2120212. 
Thus is our initial word, and Xi 1— > 2^2120212 is our substitution. 

Example 2 ((rfinBrfC) 00 ). Define a sequence as follows: * = *%* 2 . . .* n *n+\ ■ ■ ■ 
Now consider o f% o f\{*) . 

fc°fB° /!(*) = fc° /fl(*i *2 ^ *! * 2 A . . .) 

= fc(*i B *2 BAB *! 5 * 2 BAB . . .) 
= *i£ * 2 BCAB * a BC * 2 BABC . . . 

Therefore our substitution is: 
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X 2n -lX 2n ^ X 2n ^iBx 2n BC ABx 2n -\BC 'x 2n B ABC 

This substitution rule has two unknown points (x 2n ~i and x 2n ), but the seed of 
(n 2 A nB fl 4 C_)°° has only one point in it (namely A'). Thus apply (fc°fB_°fc t° 
the seed until the first 2 points are defined, as follows. 

fc°fB° fc_{A ****...) = o f R (A * A* *A* *A . . .) 

= fc(AB * BAB *B*...) 
= AB* BCAB *B... 

Thus the seed of our substitution is AB. 

3 The n operation in two dimensions. 

We have described the fl operation in one dimension. We will now extend the fl 
operation to higher dimensions and illustrate how this can be used to recreate an 
aperiodic tiling. For reasons of brevity we will merely give the definition of the 
operation in two dimensions, and a motivating example. Full details and proofs can 
be found in [REF THESIS]. 

As the reader may recall (ie figure [1]) , in the 1 dimensional case we used over- 
laying infinite strips as a physical motivation behind the D operation. Any gener- 
alization to two dimensions should be roughly equivalent to taking two unit square 
tilings on overhead transparencies, overlaying them, and reading off a new tiling 
from the result. 

Consider the following definition. 

Definition 20 (Af]B for Z 2 -arrays). Let A : Z d \-> £ and B : Z d ^ £ be Z 2 arrays 
with the standard lattice with points at integer positions. 

A n B is an array with alphabet £. The underlying lattice of the array differs 
from the standard lattice, and not in an intuitive way. 

Let A = (72' 73)' ^ 2 = 75)- 

Note that |/i| = I/2I = 1- Let LattAnB = {difi + 02/21^1,02 € Z d } be the 
underlying lattice for An B. 

Then define; 
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(AnB)(m,n) 



m—n m+n \ \ 
2 12 11 



if m + n is even 
)) if m + n is odd. 



m— re— 1 m+n— 1 
2 ' 2 



Definition 21. Let X and A be Z 2 arrays. Then define Qa{X) to be the function 
sending X to X (1 A. 

These definitions are of the same format as the definition of A D B in the one 
dimensional case. It can be rephrased as follows 

A point xe\ + ye2 in A is sent to the point xfi — yji in A fl -B. A point xei + ye2 
in B is sent to the point (x + l)/i + yfi in An .B. This covers the whole of the array 



Definition 22. Denote an array with underlying lattice generated by the vectors 
ei,e2 as a (unit) square array. Denote an array with underlying lattice generated 
by the vectors fi , / 2 as a diamond array. 

For an intuitive construction of the (AdB) operation, see the following remark; 

Remark. Take two unit square tilings, A and B such that the origin of R 2 is the 
centre of a tile in both A and B. Switch to the lattice view of tilings (ie, consider the 
array with underlying lattice formed by the centres of the unit square tiles) . Shift B 
by the vector v = (^, ^). Overlap the two lattices. The vector v has been chosen 
to maximize the distance between the set of tile centres of A, and the set of tile 
centres from B, i.e. maximizing d(A,B). From the tile viewpoint, this corresponds 
to choosing the vector v to place the centre of tiles from B at the vertices of tiles 
from A. 

Scale up so that d(B, A) — 1. The resultant array is no longer generated by the 
standard basis of R 2 . We will thus change the generating vectors, to a minimum 
vector between points in A and B, and a vector orthogonal to that minimum vector. 

This is reflected in the non-intuitive function i n ~~^~~ ' ; anc l the basis vectors e± 
and e2 changing to f\ and /2. See figure [5] for a pictorial representation of A, B and 



As the reader will be aware, An B has a different underlying lattice to A or B. 
The corresponding tiling is not a tiling by unit squares, but by unit diamonds. 

If we wish to consider multiple fl operations, it is much more natural to consider 
sequential pairs of fl operations. 



An B. 



AnB. 
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Figure 5: Patches of the tilings A, B and A D B 
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3.1 nAnB 



We already have a function (from definition I2TI) which applies an unit square 
array A to some unit square array X to produce a diamond array X fl A. 

In order to consider the limit of this qa function (and hence the Pi function), 
we must find out how to apply it to a diamond tiling. Let A and B be unit square 
arrays. Let X be some unit square array. 

Remark. Like in the ID case, we will use a null array to help clarify how the (nAf) 
B) operation builds up the periodic parts of (r\Ar\B)°° through each application of 
(C\A n B), and to show how much of the array is uniquely defined by the repeated 
application of (C\A fl B). 

In the previous section, we showed that X D A is a diamond array. Thus X fl A 
has different basis vectors to a unit square array. Let us see what effect naively 
applying the fl operation to a diamond array will have. 

From the previous chapter, the gs operation applied to any array X, will double 
the distance of every point from the origin, (and then add points from B). 
Our points in (X fl A) are of the form -^(a — b, a + b) for a, b G Z. 

Thus gs{X fl A) = 2(a — b, a + b) for a, b G Z, which, when represented as a 
tiling, is the subset shown in figure |6] (points from X and A are labelled as such). 

This is precisely half of an unit square array. Thus we need to ensure that the 
points from array B fill the gaps. If B was a diamond array, this would simply be 
a case of scaling B by V% and shifting it by a small vector. Of course, B is not a 
diamond array (though you could easily define a related fl operation where B is a 
diamond array. 

Instead we will alter how the n operation will effect points from the B array, as 
follows; 

Definition 23 (A fl B, if A is a diamond array). Let A be a diamond array, with 
basis vectors /i,/2- Let B be a unit square array, with standard euclidean basis 
vectors e\, e 2 . 



Then 




if m + to is even, 
if rro + to is odd. 



The underlying lattice for A fl B will be {n\e\ 



+ 



TO 2 e 2 |TOi, to 2 G Z}. 
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Figure 6: The image of X D A under the D operation 
This definition can be rephrased as follows. 

A point xfi + yf 2 in A is sent to the point (x — y)ei + (x + y)e 2 in A PI B. A 
point xei + ye 2 in B is sent to the point (x H — y + l)e 1 + (x + y)e 2 in A n B. This 
covers the whole of the array A n B. 

Effectively, this function takes the underlying lattice of the square array B, 
rotates it anti-clockwise by j, then moves the point of the lattice which was over 
the origin, to the point (0, 1) G Z 2 . The underlying lattice of A is then used to fill 
in the gaps in Z 2 , as in figure [71 

Obviously the rotation effect is not ideal, but its effect can be ignored by choosing 
a B tiling which has only one prototile (as we will do later). 

3.2 Motivation 

As motivation for the D operation applied to diamond tilings, let us extrapolate 
from our physical motivation with square tilings. With square tilings, we aimed to 
'overlay' the tilings, placing the second square tiling B such that every centre of a 
tile in B was as far away from the centre of tiles in A as possible. We would then 
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Figure 7: (inAflB) (schematic) 

scale up this new tiling to ensure that the minimum distance between centres of tiles 
was 1. 

The problem we will get by applying this motivation to a diamond tiling (XnA) 
and a square tiling B is that we have to rotate B to ensure the points of B are as 
far away from points of A as possible (ie, maximising d(B,A)). 

We will have several choices of how much we rotate B by, namely (? + /c|, 
k G {0, 3}), just as we have several choices of how far we translate B by (any (a, b) 
such that a + b is odd). We will choose ^ and (0, 1) in this thesis. 

We have defined X n A, when X is a diamond array and A a square array. We 
will now define C\A (1 B. 

Definition 24. Let X be a unit square tiling. Let A and B be unit square tilings. 
Then X fl A fl B is evaluated from left to right, as in the 1 dimensional case. 
Similarly longer strings of fl operations are evaluated from left to right, for example; 

x(nA n B) n = (X(nA n B) n ~ l ) n A) n B 

Also, if the limit of X(nA n B) n as n -> oo exists, denote it as X(nA n B)°°. 
We can also consider the operation as a function from the domain of unit 
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square tiling into the domain of diamond square tilings. Similar functions can be 
defined from diamond square tilings into unit square tilings (distinguishable by the 
domain specified). 

Definition 25. Define the function g% o g^{X) as fs(X R A). 

Remark. Note that HA R B doubles the distance of any point (x, y) G X from the 
origin. This is because DA sends the (x,y) to y/2(x,y), then (IB sends y/2(x,y) to 
2(x, y). Thus in the limit of (DAdB)". there will be only one tile from X remaining, 
at the origin. Thus in the current form, X R (A R B)°° is not quite Toeplitz, since 
the origin is not defined solely by f]A R B. In later examples, we will set X = B 
because of this problem, to sidestep it. 

Theorem 7. If A and B are unit square arrays, (C\A R B)°° is well defined. 

Proof. Firstly, note that X((lA R B) n is a unit square tiling, for all n. Thus all 
we need to show is that each tile t has an unique label, for all X(C\A R B) n where 
n> N t , where N t G N. 

Divide X (1 Ad B into three sets; the origin, points (x, y) where x + y is even, 
and tiles (x + y) where x + y is odd. The origin of X R A R B is unchanged under 
the operation RA R S (being the origin of X). 

The odd tiles of X R A R I? only depend on the tile types of B. Since any tiling 
X(RA R B) n can be rewritten as (X R A R B . . . R B) R A R 5, these tiles are fixed 
for all n > 1. 

Let us consider the tiles where x + y is even. The tiles in X R A R B where x 
and y are odd are uniquely defined by A. Thus, since X((]A(]B) n can be rewritten 
as (X R A R B . . . R B) R A R 5, this holds true for X(RA R £) n as well. For the tiles 
where x and y are &oi/i even, recall that the gs Qa function doubles the distance of 
every point from the origin. We know any point (x, y) where x + y is odd, is fixed 
for n > 1. Thus any point (a, b) = 2 k (x, y) where x + y is odd is fixed for n > 1 + k. 
Call this the doubling property. 

Take a point (x, y), where x and y are even. By the prime factorization theorem, 
x and y can be decomposed uniquely into factors. Thus (x, y) can be rewritten as 
(2 q s,2 r t), where s and t are odd numbers (possibly 1), and q and r are integers. 
WLOG assume q > r. Via the doubling property, we know that after r iterations 
of gs ° gAi {2 q s,2 r t) will have the same label as (2 q ~ r s,t), which is in the fixed set 
dependent on B (since 2 q ~ r s is even, t is odd, thus 2 q ~ r s + t is odd). If q — r, 
then we have (2 q s,2 r t) having the same label as (s,t), where s and t are both odd 
numbers and thus fixed by A. 
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Thus every point in X D A fl B is fixed in the limit. 



□ 



For clarification, figure M illustrates a patch of X(r\A(l B)°° , with colours added 
to indicate which tiles are defined by which iteration of qa and §b- 




Figure 8: Patch of X(nA n B) 



3.3 Illustrative example: Robinson tilings 

A Robinson tiling is a well known tiling from the field of aperiodic tilings [8] [TJ. 
A nonperiodic tiling is a tiling which is invariant only under a trivial translation. 
An aperiodic tiling is a nonperiodic tiling whose tiles cannot be rearranged into a 
periodic tiling. A sample patch of a Robinson tiling can be found in figure El 

The defining characteristic of this tiling is a hierarchy of interlocking squares. 
In this tiling, four 3x3 squares form the four corners of a larger 7x7 square. Four 
of these 7x7 squares form the corners of a 15 x 15 square, which themselves form a 
31 x 31 square, and so on. This section will show how to construct the interlocking 
square structure via a two dimensional version of the fl function. (The structure is 
'Mutually Locally Derivable' to the Robinson tiling, as per the definition in [9]). 
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Figure 9: Patch of the Robinson tiling 
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Let B be a tiling by blank unit tiles. Let A be the periodic tiling depicted in 
figure [TQl with 4 prototiles with corner decorations. Note that the tile over the 
origin has a marking with a corner in the south west of the tile. 



• marks the origin point 
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Figure 10: A 



Use B as our seed tiling. If we construct B C\ A C\ B, we produce the tiling 
depicted in figure [TTJ 

This pattern bears a certain similarity to the aperiodic Robinson tiling as shown 
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marks the origin point 



Figure 11: B(nAnB) 



in [91 More precisely, the tiles from A with corner decorations can be identified with 
cross tiles of the Robinson tiling. However only cross tiles which are corners of 3 x 3 
tiles have a counterpart in B D A fl B. 

Let us consider B(nA D B) 2 

As the reader can see, applying 9b° 9a again to B D A D B adds another layer 
of corner tiles, corresponding to the cross tiles of 7 x 7 tiles. This is due to the 
doubling property of the gs ° 9 a operation sending a 3 x 3 square onto a 7 x 7 
square. Similarly, B(nA fl B) 3 will provide corner tiles corresponding to the cross 
tiles for the 15 x 15 squares, and in general B((]A fl B) n will provide the cross tiles 
for all 2 n x 2™ squares (as well as cross tiles from smaller squares). 

Since B(jlA fl 5)°° is well-defined, in the limit you get a tiling with corner tiles 
corresponding to every cross tile in a Robinson tiling. Since the cross tiles encode 
the positions of all 2 m x 2 m squares (for m G Z) in the Robinson tiling, this is enough 
to encode the specific Robinson tiling. Thus B((lA fl 5)°° is MLD to a Robinson 
tiling. 

To be precise, it is MLD to a Robinson tiling with four faultlines meeting at the 
origin. 
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Figure 12: B(f]A n Bf 



Remark. Other aperiodic tilings similar to the above Robinson tiling can be gen- 
erated by changing the origin point of gs ° Qa- For further details, see [3]. 

Remark. The two dimensional variation of the D operation bears similarities to 
certain areas of research on square tilings. In particular this concept may have uses 
with the Winfree model [H] (self-assemblying constructions based on square tiles), 
since we have seen that the method can construct at least one two-dimensional 
structure. 

Another possible use for the D operation is in the field of computer graphics. 
Methods to create computer graphics from Wang tilings are well studied and form 
one of the most common methods of texture generation (for example [Uj). Briefly, 
a simple texture is placed onto each Wang prototile, and a tiling is formed from 
the Wang prototiles. The eye cannot see any periodicity in the tiling, thus a con- 
vincing graphic can be produced from simple textures. This method breaks down 
when the graphic is approached, since the original simple textures can be seen more 
clearly, and pixelation may occur. Methods to stop this problem are currently being 
investigated (for example [ID]). 

The two dimensional D function lends itself to this problem, since if you remove 
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the rescaling factor of the D function, the size of the tiles shrinks with every iteration. 
Thus by applying the D function multiple times as the graphic is approached, new 
tiles can be introduced, and thus pixelation should be avoided. 

On a more mathematical note, the last two chapters have described operations 
from limits of D operations to Toeplitz sequences. Investigating whether a reverse 
operation can be constructed from a subset of Toeplitz sequences to compositions 
of H operations would be an interesting avenue of inquiry. 
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